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Abstract. The interaction of matter- wave solitons with a potential barrier 
is a fundamentally important problem, and the splitting and subsequent 
recombination of the soliton by the barrier is the essence of soliton matter- 
wave interferometry. We demonstrate the three-dimensional (3D) character 
of the interactions in the case relevant to ongoing experiments, where the 
number of atoms in the soliton is sufficiently large that the soliton is close 
to collapse. The mean-field description is quite accurate, but the proximity 
to the collapse threshold makes the use of the ID Gross-Pitaevskii equation 
(GPE) irrelevant. We examine the soliton dynamics in the framework of the 
effectively ID nonpolynomial Schrodinger equation (NPSE), which admits the 
collapse in a modified form, and in parallel we use the full 3D GPE. Both 
approaches produce results which are very different from those produced in recent 
work which used the ID cubic GPE. The features observed in both the NPSE 
and the 3D GPE include (a) an increase in the first reflection coefficient for 
increasing barrier height and decreasing atom number; (b) strong modulation 
in the secondary reflection/recombination probability vs. barrier height; (c) a 
pronounced asymmetry in the oscillation amplitudes of the transmitted and 
reflected fragments; (d) an enhancement of the transverse excitations as the 
number of atoms is increased. We also explore effects produced by variations of 
the barrier width and outcomes of the secondary collision upon phase imprinting 
on the fragment in one arm of the interferometer. 
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1. Introduction 

Ongoing studies of atomic Bose-Einstein condensates (BECs) have contributed 
numerous fundamental insights in a wide range of phenomena [TJ [2]. This may 
be largely attributed to the precise control afforded by experiments, and to the 
existence of accurate, yet quite tractable models based on the Gross-Pitaevskii 
equation (GPE) [3JH]. Many of the theoretical and experimental investigations have 
strong connections to other areas, such as condensed matter physics, nonlinear optics, 
superconductivity, and superfluidity. Furthermore, ultracold gases trapped in external 
potentials may be utilized as quantum simulators of real materials [5.J. 

Inter-atomic interactions in a BEC enable the examination of a wide variety 
of nonlinear effects that have been summarized in a number of reviews and books 
such as [4]. Basic types of coherent structures built by the nonlinearity are bright 
solitons and soliton complexes [6j [TJ El El HH] , dark solitons [TT] . vortices, and vortex 
lattices [T2l [13j HH [15] . Bright solitons were proposed for matter- wave interferometry 
[16], which is itself an active field of research [17 . Recent work focused on the role 
of the condensate's effective nonlinearity [18 , as well as on advanced applications 
[T9] . Solitons may provide for 100-fold improved sensitivity for atom interferometry, 
where their long lifetime, of the order of seconds, may enable precise force sensing [20] 
and measurement of small magnetic field gradients [21 . Robust bright solitons, 
appropriate for these applications, have been created under well-controlled conditions 
in BECs with attractive [TBI [22j [23] and repulsive [24] inter-atomic interactions (in 
the latter case, these were gap solitons [7] supported by the optical-lattice potential). 
A topic of particular interest has been the splitting and subsequent recombination 
of solitons due to collisions with potential barriers. This was studied theoretically 
in detail [25j [26j [271 EH], as summarized in a recent review [29 . However, it is 
relevant to stress that the analysis has been thus far carried out with a one-dimensional 
(ID) model, i.e., the ID GPE with cubic nonlinearity. In the same framework, the 
quasiparticle-wave duality of solitons and their tunneling through a potential barrier 
was recently considered [31] [32] . 

The subject of the present work is the interaction of the incident soliton with 
a potential barrier, acting as a beam splitter, and the subsequent recombination of 
the split fragments in the 3D setting. The three-dimensionality strongly affects the 
results in comparison with the ID setting of the above recent studies. Bright matter- 
wave solitons are experimentally created in cigar-shaped (nearly-lD) potential traps 
in which the self-trapping is driven by the intrinsic attractive nonlinearity solely in the 
axial direction. The number of atoms in the soliton, N will be considered close to the 
value at the collapse threshold, 7V C , and the actual shape of the solitons will be three- 
dimensional in the considerations that follow. Accordingly, the solitary waves do not 
feature the conspicuous elongation of the underlying trapping potential [I6j [22J [23] . 
Consequently, the character of the splitting and subsequent merger may be essentially 
three-dimensional — in particular, transverse modes may (and will, as we show) be 
excited. 

In the recently studied ID counterpart of this setting, results of an earlier work [33] 
were used in [26 to predict the transmission and reflection coefficients and the phase 
shift resulting from the interaction with a narrow barrier. A sinusoidal dependence of 
the transmission coefficient on the relative phase of two solitary waves colliding on the 
barrier was identified. In [27J, both the primary splitting of a single solitary wave and 
the secondary collision of the two fragments was studied using the truncated Wigner 
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approximation, with emphasis on the comparison of the mean-field dynamics with 
the quantum dynamics. Significant deviations between the two were identified for a 
relatively small number of atoms in the solitary wave. A step-wise (discontinuous) 
variation of the reflection and transmission coefficients due to quantum superposition 
was recently reported [28j [30j [31] . The deviation between quantum and mean-field 
approaches is significant only for relatively small TV, becoming negligible for large TV. 
In this paper, we consider TV close to collapse, where TV C ~ 10 5 . For this reason, 
quantum effects are not crucially important to us. 

Collisions of ID solitons with an attractive potential well, rather than the 
repulsive barrier, have also been studied recently, revealing a fairly complex 
phenomenology. The latter features alternating windows of transmission, reflection, 
and trapping with abrupt transitions between them as the potential depth varies [3H 
[35]. 

Our analysis is performed in two stages. First, we use the well-known ID 
nonpolynomial Schrodinger equation (NPSE) [36l [37], which, unlike the cubic ID 
GPE, captures the 3D phenomenology in some approximation (in particular, it admits 
the possibility of collapse, although not in the same way as the 3D GPE). We will 
refer to this NPSE setting as a quasi- ID description to clearly distinguish it from 
the ID GPE setting where collapse is absent. Secondly, we report results obtained 
from the full 3D GPE. To the best of our knowledge, both equations are used for 
the first time in this context. Our main interest is to consider relatively large atom 
numbers TV (TV = 80000 ~ 0.85TV C ), for which the phase coherence of the condensate is 
maintained and the soliton dynamics are most interesting. We explore the dependence 
of the results on TV, as well as on the height of the barrier, E. Between the dynamics 
described by the ordinary cubic ID GPE and predictions of the full 3D equation, there 
exist significant differences due to the fact that of exciting transverse oscillations (both 
absent in the ID description), features that become more prominent as TV — >> TV C . On 
the other hand, the quasi-lD NPSE is found to be in reasonable agreement with the 
3D GPE both in regard to its static properties (although the stable branch of solutions 
at low TV is captured more accurately than the unstable collapsing one of high TV) and 
to the dynamics of soliton-barrier interactions. 

The paper is structured as follows. In Section 2, we present the quasi-lD and 
3D models used in the computations herein, compute their stationary solutions, and 
define the diagnostics used to examine the dynamics. In Section 3, we report the 
numerical results in detail, varying the barrier height and width, atom number, as 
well as phase imprinting applied to one of the fragments to simulate interferometry. 
The paper is concluded by Section 4. 



2. Model equations and stationary solutions 

We start by introducing the quasi- ID NPSE for the mean-field wave function, ip, of 
the condensate of 7 Li atoms loaded into a cigar-shaped trap with axial coordinate z, 

dip h 2 d 2 ^ i + (3/2)gM 2 

in—- = — - — — — + V{z)w H —ib, (1) 

dt 2m dz 2 v ^ yrr^w 

where g = 2a s /a_\_, with the typical scattering length a s = —0.3 ao, where ao is 
the Bohr radius, and transverse trapping radius a_\_ = y/h/(mu)±) = 2.25 /im, which 
corresponds to a confinement frequency uj± = 27T-290 Hz and atomic mass m = 7 amu. 
The longitudinal potential V(z) includes a weak parabolic trap, which we assume has 
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a typical value of the strength, uj z = 2n • 5.6 Hz, and the Gaussian barrier of height E 
and width e, which will be varied below: 

V(z) = (l/2)muj 2 z z 2 + £exp (-2z 2 /e 2 ) . (2) 

The total number of atoms, given by N = ^ IV ; ( 2: )| 2 ? wml a l so be varied. 

The results produced by the NPSE u\l will be compared to those obtained from 
the radially symmetric 3D GPE, written in the cylindrical coordinates, (p, z) as 

M h 2 (d 2 ^ l top d 2 i)\ inn 2 2 



with N = 2tt J °° pdp dz \ip (p, z)| 2 and the 3D potential, 

V(p,z) = (l/2)m(uj ± p 2 +uj 2 z z 2 ) + £exp (-2z 2 /^ 2 ) • (4) 

Equation Q implies that we consider only fundamental axially symmetric solitons, 
but not ones with intrinsic axial vorticity, which are also possible in this setting [37 . 
Axial widths of the solitons generated by the quasi- ID NPSE and 3D GPE are defined 

as 



/ r+oc / poo />+oo 

W 1D = 2JN- 1 j dzz 2 \t/>\ 2 , W m = 2 JN- 1 J pdp J dzz^. (5) 

Stationary solutions with chemical potential p were sought for as ip(r,t) = 
ex.p(—ipt/h)ijj(r). Figure [l] compares N(p) and W(N) curves for the soliton families, 
produced by both the quasi- ID and 3D equations. In these plots, N is normalized to 
the critical atom number, N C1 which initializes the collapse, i.e., the largest number 
atoms accessible (by a stable standing wave solution) in the given setting [38j [39] . 
Both equations yield nearly the same value, 7V C w 94000. 

It is observed that stable solution branches generated by the NPSE and GPE 
are very close, while the discrepancy is conspicuous for the unstable branches. In the 
right panel of Figure [l] a horizontal line is drawn at W = 4.5 pm, which corresponds 
to the typical width of the barrier e in the simulations reported below. It is seen that 
stable solitons are narrower than 4.5 pm only at N > 0.987V C . Consequently, solitons 
are typically wider than the considered barrier. 



3. Numerical simulations 



3.1. Effect of changing the number of atoms 

We first consider simulations for the solitons launched with zero velocity at the distance 
A = 251 pm from the central position, where the barrier is placed. In the present 
setting, this amounts to the initial potential energy equal to 678 Hz. The barrier is 
here taken with a fixed width, e = 4.5 pm. The number of atoms in the soliton and 
the barrier height E are varied and the outcome of the soliton-barrier interaction is 
monitored. The splitting is characterized by the time-dependent reflection coefficient, 

/0 rOO rO 

dz\x[>\ 2 , R 3D (t) = N~ 1 pdp dz\ip\ 2 , (6) 
-oo •' J — oo 

which shows the fraction of atoms remaining in the incident arm of the interferometer. 
Below, we focus on R taken at two specific times, R\ = R(t = tt/uj z ) and R2 = R(t = 
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Figure 1. (a) The number of atoms, normalized to its critical value at the 
collapse threshold, N/N c , versus the chemical potential, for the stable (/i > /i c ) 
and unstable (/i < n c ) soliton families; /i c = 147.9 Hz in the quasi-lD NPSE 
setting |l} and fi c = 130.5 Hz within the 3D GPE (3}. (b) The axial widths, 
denned as per j5}, versus N/N c , for the stable soliton family. The green dashed 
horizontal curve of the right panel corresponds to the width of 4.5 /im, typically 
associated with the barrier in what follows. 



2tt/uj z ). The former value measures the fraction of atoms that remain in the original 
arm after the interaction of the incoming soliton with the barrier, because t = tt/uj z 
corresponds to a half-period of oscillation of the soliton in the weak longitudinal trap, 
thus ensuring that its first collision with the barrier (occurring roughly at a quarter 
of the period) is completed. The latter value, measures the relevant fraction after 
the second interaction, i.e., the collision of the returning fragments with the potential 
barrier. This occurs at roughly three quarters of the oscillation period, while R2 is 
measured at a full period (defined in the absence of the barrier). 

Figures [2] and [3] show, respectively, R\ and R2 versus the number of atoms, TV, 
and the barrier height, E, with the blank regions denoting the presence of collapse. 
Collapse occurs for N < N c in these cases due to the interaction with the barrier. It is 
observed that the quasi- ID NPSE and 3D GPE dynamics show similar trends for the 
first reflection. Namely, the reflection increases with the growth of the barrier height 
and with decreasing N. The second reflection coefficient, i?2, presents a more complex 
functional dependence, which, in the terms of the quasi-lD NPSE and 3D GPE alike, 
features an oscillatory variation with N and E for sufficiently large N . Notice the very 
good qualitative and even semi-quantitative agreement between the quasi- ID and the 
fully 3D models in this context. This behavior can be attributed to the number- 
dependent variation of phases and amplitudes of the fragments emerging from the 
first collision, which, upon their recombination, leads to outcomes ranging from the 
nearly perfect reflection of the recombined soliton to its nearly perfect transmission. 
This is a more complex manifestation of the feature observed in Ref. 26] and is a 
nonlinear effect stemming from the interference of the phases of the two fragments, 
becoming progressively more pronounced as N — >> N c . This phenomenon is especially 
visible for N/N c > 0.3, leading to a series of resonant peaks of full reflection. In 
some sense, this effect may be considered as a nonlinear analog of the Fano resonance, 
which is well known in linear wave systems [40]. Another interesting feature is the 
presence of finger-like gaps in the quasi- ID setting where collapse occurs after the 
second interaction with the barrier, when N/N c > 0.85. Those gaps do not appear 
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Figure 2. The first reflection coefficient, R\ [see |6f] vs. barrier height, E, and 
the normalized number of atoms, N/N c , in the quasi-lD (a) and 3D (b) settings. 
Blank regions correspond to barrier-induced collapse. Panel (c) shows the value 
of the barrier height at which equal splitting (Ri = 0.5) is attained, as a function 
of N/N c . The barrier width is e = 4.5 /im. 



in the 3D setting, indicating that they originate from an inherent limitation of the 
NPSE approximation. This is explained by the fact that in the framework of the 
NPSE the soliton collapses when the denominator vanishes, at the critical value 
of the peak density (contrary to the genuine collapse of the fully 3D setting). This 
difference also justifies the apparently wider range of parameters (in the (E,N/N C ) 
parametric plane) leading to collapse for the quasi- ID problem in comparison to the 
fully 3D one. Finally, the bottom panel of Figure [3] shows the increase in the second 
reflection with increasing iV, for the case in which E is adjusted to give 50% splitting 
after the first reflection. 

While the reflectivities provide a measure of the asymmetries between the atom- 
number fractions emerging to the left and right of the barrier, we have used additional 
diagnostics to quantify the dynamics. In particular as a measure of the asymmetry of 
the reflected and transmitted waveforms, we define 

A t +Ar' 

where A r and A t are the oscillation amplitudes within the trap of the reflected and 
transmitted fragments, and are both taken to be positive. The dependence of f on 
N/N c and E is presented in Figure [4j along with its dependence on N/N c when E is 
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Figure 3. Panels (a) and (b) depict R2 vs. E and N/N c , otherwise the same as 
Figure [2] The blank areas indicate where barrier-induced collapse arises. Panel (c) 
shows the typically increasing trend of R2 as a function of N/N c when R\ = 0.5. 
The barrier width is e = 4.5 /im. 



adjusted to give 50% splitting. An interesting feature is that £ > for every N and 
E, i.e. the transmitted fragment always reaches a higher \z\ and, consequently, has a 
higher kinetic energy, independently of whether it is the larger or smaller fragment. 
This can presumably be attributed to the original direction of motion (and associated 
momentum) of the incoming solitary wave. It should also be noted that typically 
the most pronounced asymmetries occur for larger values of N/N c , appearing to be 
predominantly a feature of the 3D nature of the interactions in that case. 

To assess the impact of the width of the barrier on the above results, we have 
performed a complementary study for the ID NPSE, considering a narrower barrier 
whose width is e = 1 /im, a value much less than the axial width of the soliton. Figure 
[5] shows the dependence of R\ (a), R2 (b) and ( (c) with respect to N and E as before. 
The qualitative nature of the dependence of these quantities does not seem to change 
in comparison to the case of e = 4.5 jam. Nevertheless, there is a quantitative shift of 
the principal features towards higher values of the barrier height E, as expected. 

A basic feature of the 3D dynamics, which may be naturally expected, and is 
indeed produced by the 3D GPE, is the excitation of radial oscillations after splitting. 
This feature is, by construction, not incorporated in the ID GPE equation, where the 
transverse shape is assumed to be "frozen" . On the other hand, it is incorporated in a 
simple (yet reasonably accurate) way in the quasi- ID description of the NPSE, through 
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Figure 4. (a) Oscillation asymmetry, £ [see the definition of j7|], after the 
first barrier interaction as a function of E and N/N c , in the 3D setting (as with 
previous panels, the contour plot is essentially identical for the quasi-lD case of 
the NPSE). (b) C versus N/N c when Ri = 0.5. The barrier width is e = 4.5 /im. 
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Figure 5. First reflection coefficient, R\ (a), second reflection coefficient R2 (b), 
and the oscillation asymmetry C, (c) vs. E and N/N c , in the quasi-lD setting of 
the NPSE when e = 1 /im 
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Figure 6. (a) Amplitude po of radial oscillations in the 3D setting excited by 
the first interaction with the barrier, as a function of E and N/N c . In (b), the 
value of po is plotted versus N/N c , while R\ = 0.5 is fixed. The barrier width is 
e = 4.5 jum. 



the assumption of a space- and time-dependent width of the transverse direction of the 
ground state which is directly controlled (according to the Euler-Lagrange equations) 
by the longitudinal wavefunction. This feature is partially responsible for the good 
agreement between the quasi- ID NPSE and the fully 3D GPE, as observed above. 
We quantify the radial vibrations in the 3D case by evaluating a measure of their 
amplitude 

m = Jf Jo pdp J dzp ^ 2 > ^ 

and defining po as the maximum value of p (over time for a given set of parameters). 
Figure [6ja) shows the dependence of po on N/N c and E, together with its value at 
the barrier height E sp corresponding to R\ = 0.5 (b). It is clearly observed that with 
the increase of N/N c and E, the excitation of transverse oscillations becomes very 
significant, with the amplitude attaining values ~ 0.5 pm. The same conclusions are 
suggested by the right panel for the case of even splitting. Hence, it can be inferred that 
the role of higher-dimensionality (captured qualitatively within our quasi- ID NPSE 
approach and properly incorporated in the fully 3D setting) is of particular relevance 
for iV's close to N c (a case of central interest to ongoing experiments). 

Figures [7] and [8] show outcomes of simulations corresponding, respectively, to 
small and large values of i?2, in the 3D setting (results for the quasi-lD NPSE are 
quantitatively similar). The figures display the longitudinal (z-dependent) density, 
resulting from the integration of the density in the transverse plane (denoted as n z ), 
along with amplitudes of the fragments and the time-dependent reflection coefficients, 
which exhibit jumps upon the interaction of the one (during the first event) or of the 
two (during the second event) solitary wave(s) with the barrier. As these figures show, 
the amplitudes of the fragments can be different even when the numbers of atoms in 
them are equal. Naturally, the width is larger for the fragment with the lower peak 
density. Furthermore, in addition to the two limiting cases, our results suggest that 
one can manipulate the parameters (such as E, A/", etc.) controlling the interaction 
of the incident soliton with the barrier to produce any desired outcome within a wide 
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Figure 7. Interaction with a barrier of width e = 4.5 /im giving Ri =0.5 and 
small R2 (i.e., even splitting of the incident soliton followed by the recombination 
of the fragments into a nearly single transmitted one) in the 3D GPE (a) 
Integrated density plots, (b) The amplitude of the fragments divided by the 
amplitude of the initial soliton (c) The time-dependent reflection coefficient. The 
present results were obtained with N « 0.17iV c and E = 696.95 Hz. The smallest 
value of the second reflection coefficient is R2 = 0.17. The outcome for the quasi- 
1D NPSE is essentially similar. 



range in the sequence of two collisions. While there appear to be two waves in Figure 
[7] for z < 0, after the second soliton-barrier interaction, this feature can be controlled 
(and avoided) through the use of a narrower barrier (results not shown here). 

3.2. Effect of changing the initial position of the soliton 

In this subsection the number of atoms is fixed to N = 0.87V C , which is a typical, 
experimentally reachable value where the 3D nature of the condensate is crucial to 
the observed dynamics. Hence, we expect the 3D GPE to be a suitable model, and the 
quasi- ID NPSE to be a good qualitative (and reasonable quantitative) approximation. 
We also expect transverse excitations to be present in the dynamics. 

We display here the dependence on initial displacement A from which the soliton 
is launched with zero velocity. For completeness, we do this for different barrier widths 
m (Q, e = 1,2,3,4 and 5 /im, in the 3D setting. Figure [9] shows the dependence on A 
for the barrier height E sp giving even splitting, the second reflection coefficient (R2), 
the asymmetry parameter ([7]), and the radial oscillation amplitude, po. As expected, 
Figure [9la) shows that E sp increases with increasing A. Under these conditions of 
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Figure 8. The same as in Figure [7] but for the case of R2 close to 1, i.e., the 
recombined soliton returning in the incoming arm. iV « 0.887V C and E = 686 Hz, 
resulting in Ri = 0.44 and R 2 = 0.96. 



strong nonlinearity (N = 0.87V C ), however, E sp depends only weakly on the barrier 
width £, with the exception of the narrowest barrier. Figure [9jb) shows that R2 
has the least sensitivity to A for the narrowest barrier with the required height being 
larger for smaller width, while R2 is larger (as may be expected intuitively) for a larger 
width. An increase in A generally leads to only a weak modification of the asymmetry 
factor £, which is chiefly decreasing with A for wider barriers, presumably due to the 
larger speed of the soliton impinging upon the barrier, and correspondingly smaller 
interaction times. For narrow barriers, the dependence on A is more pronounced 
and non-monotonic. Figure [9^d) shows that transverse excitation is stronger for more 
energetic solitons and narrower barriers, where the interaction with the barrier is more 
impulsive. We expect transverse excitation to be more probable when the atomic 
kinetic energy is comparable to the transverse mode spacing huoj_. This occurs when 
A = 164/im, which is reasonably consistent with the sudden increase in p observed 
for A ~ 100/im. 

We note that the most significant, although small, differences between the quasi- 
1D setting of the NPSE and the full 3D setting arise in the case of narrow barriers. 
This is intuitively reasonable, as a narrow barrier induces dynamics on length scales 
closer to the transverse confinement and hence, enhances the degree of transverse 
excitations thus affecting the quality of the approximation of the 3D behavior by the 
quasi- ID NPSE. Nevertheless, the agreement between the two is still in reasonable 
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Figure 9. Dependencies of characteristics of the interaction on the initial 
displacement of the incident soliton, A, for N = 0.8N C . (a) The barrier height, 
E, necessary for 50% splitting after the first collision. For the following plots, 
E = E sp : (b) The secondary reflection coefficient, R2] (c) The asymmetry 
coefficient C, defined as per jT]); and (d) Amplitude of the radial oscillations, 
po- The displayed results are produced by the 3D GPE (3}. 



qualitative agreement. 

3. 3. Effects of the phase imprinting 

Finally, we briefly analyze the effect of imprinting a phase difference at a certain time, 
t = tt/uj z onto the initially transmitted fragment, thereby emulating the operation 
of an interferometer. The phase difference is introduced when the fragments are 
located at the largest distance from the barrier. Figure [10| shows the second reflection 
coefficient, i?2, in the cases of small (for comparison) and large atom numbers, as 
obtained from the simulations of the quasi-lD NPSE ^ and the 3D GPE It 
is seen that, in the regimes of small TV, the absence of a phase shift provides for a 
nearly complete transmission, while a relative phase shift of ir leads to almost perfect 
reflection. The contrast is considerably reduced for strong nonlinearity and sufficiently 



wide barriers, as shown in the right panel of Figure 10, although a wider range of 
transmittivities/reflectivities is accessible for narrower barriers. In that light, although 
highly nonlinear waves (such as solitons) may be deemed less useful for interferometric 
purposes, their relevance may be (at least partially) restored in the context of suitably 
narrow barriers. 
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Figure 10. Second reflection coefficient, R2, in the case of the initial 50% 
splitting, as a function of the phase difference introduced between the fragments 
at the time of their largest separation. The panels correspond to (a) N « 0.067V C 
(near the linear limit, shown for comparison) and (b) TV « 0.87V C . Notice that in 
the former case, the results for ID and 3D settings coincide for e = 4.5 /im. In 
all cases, the initial displacement of the soliton was A = 251 fim 



4. Conclusion 

We have explored the collision of matter-wave solitons with a central barrier of 
variable width and height inserted into a shallow harmonic-oscillator axial trap, and 
the subsequent collision of the split fragments upon their return to the defect. This 
configuration forms the basis of a soliton interferometer. The non-monotonic variation 
of R2 with N and E, evident for strong nonlinearity in Figures [3] and [5j is detrimental 
to the sensitivity of an atom interferometer. Figure [lOj however, indicates that the 
sensitivity to phase variation can be regained by smaller nonlinearity and that the 
adverse effect of strong nonlinearity can be partially mitigated by a narrow barrier. 
Narrow in this context, evidently means in comparison to the axial size of the soliton. 

Our analysis accounted for relatively large atom numbers, which are not much 
smaller than the collapse threshold. This strong nonlinearity necessitates a full 3D 
solution, in contrast to the ID GPE setting studied in previous works. The analysis 
was carried out, in parallel, in the framework of the quasi-lD NPSE and the full 3D 
GPE, indicating a good qualitative and even reasonable quantitative (at least not too 
close to TV = N c ) agreement between the two. A detailed computational map of the 
ensuing phenomenology has been generated as a function of number of atoms in the 
soliton, its initial distance from the barrier (which determines the collision velocity), 
and the height and width of the splitting barrier. Additionally, the effect of a phase 
shift imposed on the fragments at the moment of the largest separation was also 
examined. 

A number of general conclusions, obtained in the framework of the quasi- ID and 
full 3D settings, and their similarities and differences have been reported. While the 
results are similar between these two cases, they are essentially different from those 
generated by the ID cubic GPE, which was used previously. The quasi- ID NPSE and 
3D GPE produce similar values for the first reflectivities, systematically increasing 
with increasing barrier height E and weakly decreasing with increasing number of 
atoms N . The second reflectivity (corresponding to the collision of the fragments after 
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the first interaction with the barrier) oscillates strongly from complete reflection to 
high transmission as function of barrier height (and atom number). These oscillations 
are a fundamentally nonlinear effect, most pronounced for large N. The excitation 
of the transverse breathing mode was also characterized in the framework of the 3D 
GPE. This effect was found to become progressively more significant as the critical 
number of atoms, corresponding to the onset of the collapse, was approached. Phase 
imprinting was found to play a critical role in the outcome of the second collision, 
especially for small nonlinearity, modifying it between nearly complete transmission 
and full reflection. The highly nonlinear realm seems less sensitive to such variations, 
but the sensitivity may be restored for narrower barriers. 

It is important to corroborate these findings experimentally. Regarding further 
theoretical analysis, a challenging problem is to study deviations of the results from 
the mean-field approximation in the 3D geometry as a function of the atom number 
TV, in analogy to the recently emergent studies in the ID setting [27[ [28l [30]. An 
additional issue of interest is the use of potential wells, rather than barriers, for which 
more complex phenomenology may be expected in the 3D setting [34j [35] . Lastly, 
extending such considerations to dark solitons in condensates with the self-repulsive 
nonlinearity [41] , and to multi-component systems [32j [42] are also compelling topics 
for further investigation. 
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